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ABSTRACT 

In this paper we develop a theory of disturbances induced by the stellar tidal 
field in a fully convective slowly rotating planet orbiting on a highly eccentric 
orbit around a central star. In that case it is appropriate to treat the tidal 
influence as a succession of impulsive tidal interactions occurring at periastron 
passage. For a fully convective planet mainly the I = 2 fundamental mode of 
oscillation is excited. We show that there are two contributions to the mode 
energy and angular momentum gain due to impulsive tidal interaction: a) 'the 
quasi-static' contribution which requires dissipative processes operating in the 
planet; b) the dynamical contribution associated with excitation of modes of 
oscillation. These contributions are obtained self-consistently from a single 
set of the governing equations. We calculate a critical 'equilibrium' value of 
angular velocity of the planet £l cr it determined by the condition that action of 
the dynamical tides does not alter the angular velocity at that rotation rate. 
We show that this can be much larger than the corresponding rate associated 
with quasi-static tides and that at this angular velocity, the rate of energy 
exchange is minimised. We also investigate the conditions for the stochastic 
increase in oscillation energy that may occur if many periastron passages are 
considered and dissipation is not important. We provide a simple criterion for 
this instability to occur. 

Finally, we make some simple estimates of time scale of evolution of the 
orbital semi-major axis and circularization of initially eccentric orbit due to 
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tides, using a realistic model of the planet and its cooling history, for orbits 
with periods after circularization typical of those observed for extra-solar plan- 
ets P bs <J Sdays. Quasi-static tides are found to be ineffective for semi-major 
axes > 0.1 Au. On the other hand, dynamic tides could have produced a very 
large decrease of the semi-major axis of a planet with mass of the order of the 
Jupiter mass Mj and final period P b s ~ 1 — 4.5days on a time-scale < a few 
Gyrs. In that case the original semi-major axis may be as large as > 10 2 Au. 

For planets with masses > 5Mj dynamical tides excited in the star appear 
to be more important than the tides excited in the planet. They may also, in 
principle, result in orbital evolution in a time less than or comparable to the 
life time of the planetary systems. Finally we point out that there are several 
issues in the context of the scenario of the circularization of the orbit solely 
due to dynamic tides that remain to be resolved. Their possible resolution is 
discussed. 

Key words: tides, binary star, stellar rotation, planetary systems: formation 

1 INTRODUCTION 

The recently discovered 100 or so extra-solar giant planets have masses in the range 0.12 — 11 
Jupiter masses. They may be found at distances of several Au or close to the central star 
with periods of a few days. High orbital eccentricities are common (e.g. Mayor & Queloz 
1995; Marcy & Butler 1998; Marcy et al 2000). 

It has been suggested that giant planets may form through gravitational instability in 
a disc at large radii (e.g.. Boss 2002) or through the 'critical core mass' model in which a 
critical ~ 15 M ffi core is formed in a disc by accumulation of solids and then undergoes rapid 
gas accretion (see e.g. Papaloizou, Terquem & Nelson 1999, for a review and appropriate 
references). In this case it is expected that the cores of gas giant planets should begin to form 
beyond a radius of r ~ 4 Au, the so-called 'ice condensation radius' where the existence of 
ices facilitates the accumulation of solids. 

In order to explain the existence of the closely orbiting extra-solar giant planets one is 
then led to propose orbital migration. In principle this may occur through the gravitational 
interaction between a protoplanet and the protostellar disc (e.g. Lin & Papaloizou 1986, 
Nelson et al 2000 and references therein) or through mutual gravitational interactions among 
a strongly interacting system of protoplanets, and further circularization by tides. In this 
Paper we explore the second possibility. 
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The presence of high orbital eccentricities amongst extra solar planets is suggestive of a 
strong orbital relaxation or scattering process. The gaseous environment of a disc may act 
to inhibit such interactions until it is removed. Gas free dynamical interactions of coplanar 
protoplanets formed on neighbouring circular orbits have been considered by Weidenschilling 
& Mazari (1996) and Rasio & Ford (1996). These may produce close scattering and high 
eccentricities but the observed distribution of extra solar planets is not reproduced. 

Papaloizou & Terquem (2001) investigate a scenario in which 5 < N < 100 planetary 
objects in the range of several Jupiter masses are assumed to form rapidly through frag- 
mentation or gravitational instability occurring a disc or protostellar envelope on a scale of 
Rmax = 100 Au. If these objects are put down in circular orbits about a solar mass star, at 
random in a volume contained within a spherical shell with inner and outer radii of 0AR max 
and R m ax respectively, the system relax on a time-scale ~ 100 orbits and the final state of 
the system is independent of details of initial conditions. In fact, the evolution is similar to 
that of a stellar cluster, most objects escape leaving at most 3 bound planets. Close encoun- 
ters or collisions with the central star occurred for about 10% of cases in which a planet 
goes through a phase where it is in a highly eccentric orbit that has close approaches to the 
central star. 

These simulations naturally lead us to a situation where tidal interaction leading to 
orbital circularization and to the formation of a very closely orbiting giant planet may 
become a possibility. In this Paper we analyse the tidal interactions in detail. For simplicity 
we neglect the influence of gravitational interactions between the planets and consider only 
the tidal interactions between the planet and the central star. 

We develop the theory of disturbances induced by tides in a giant planet in Section 2. We 
consider the case of a highly eccentric orbit and a fully convective slowly rotating planet tak- 
ing into account all the specific features of this situation. This is such that the theory of tidal 
disturbances is sufficiently simple that it is possible to obtain compact analytical expressions 
for the energy and angular momentum gained by the planet after a close encounter with the 
star, see equations (44,45,47) for quasi-static tides associated with frictional processes in the 
planet and (60-65) for dynamic tides associated with excitation of the oscillation modes. 

It has been mentioned by a number of authors (Press & Teukolsky 1977, hereafter PT, 
Kochanek 1992, Mardling 1995 a,b) that when there are many periastron passages, energy 
can flow backwards and forwards between the orbit and oscillation mode. This happens in a 
situation where the dissipation time of the normal mode is larger than the orbital period of 
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the planet. The direction of energy transfer resulting from a particular periastron passage 
depends on the phase of the pulsation at the moment of periastron passage. We address 
this multiple passage problem in Section 2.9 and show that when the semi-major axis of 
the orbit is sufficiently large, the energy transfer is highly irregular as has been previously 
found numerically (Mardling 1995 a,b). Here we find a simple semi-analytic criterion for 
such irregular behaviour. When this criterion is fulfilled the mode energy grows on average 
proportional to the number of periastron passages and we can use calculations based on 
'naive' adding of the energy gained from the orbit by the normal mode in order to make 
estimates of the evolution time-scale. The similar situation occurs if one simply assumes that 
the energy transferred to normal modes is dissipated between successive periastron passages. 

In Section 3 we adopt the above approach to calculate time-scales required for significant 
orbital circularization due to tides. Finally in Section 4 we discuss our results and present 
conclusions. 

2 THE TIDAL INTERACTION OF A ROTATING CONVECTIVE 
PLANET IN A WEAKLY BOUND ORBIT 

In this Section we consider the tidal perturbations induced in a rotating planet after a 
parabolic encounter with a central star. We derive expressions for the energy and angular 
momentum gained by the planet as a result of such an encounter. By considering successive 
encounters, the results can be used to discuss planets on highly eccentric orbits as well as 
parabolic ones. In the absence of explicit dissipation only the excitation of normal modes 
or dynamic tides play a role. For a non-rotating planet the energy gain has been calculated 
by PT. The equations for the perturbations induced in a rotating object have also been 
discussed (e.g. Lai 1997 and references therein). However, only dynamic tides were discussed 
and the treatment of viscous terms in the equations of motion was incorrect. Stationary 
or quasi-static tides for very eccentric orbits have been considered by Alexander (1973), 
Hut (1981), Zahn (e.g. Zahn 1989 and references therein) and more recently by Eggleton, 
Kiseleva and Hut (1998) using a very different approach based on calculating the phase lag 
of an assumed quasi-static tidal bulge. Here we demonstrate how the expressions for energy 
and angular momentum gain associated with the the quasi static and resonant dynamic 
tides follow from a single formulation of the problem leading to one governing equation. 
In general, the linearised equations governing small perturbations of a rotating object 
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are non-separable. However, for a fully convective planet and sufficiently slow rotation the 
angular velocity Q r is much smaller than the internal 'dynamical' frequency fi* = y^^i 
where m p i is the mass of the planet. Then it is possible to construct a self-consistent theory 
of perturbations based on an expansion in the small parameter f2 r /fi* *. Here, only zeroth 
and first order terms in that expansion will be taken into account. 
We introduce the inner product (e.g. PT) 

<€il£ 2 >= Jd 3 xp(^.^), (1) 

where (£ x • £ 2 ) is the scalar product of two vectors, and p is the density. For simplicity we 
consider only the / = 2 component of the tidal perturbation which dominates for large orbital 
separations. Further since g modes are absent in the fully convective objects we consider and 
the frequencies of the p modes are always significantly larger than the frequency of the / 
mode making their tidal excitation relatively small, we take account only of the fundamental 
mode when calculating the tidal response, t. 

2.1 Basic equations 

We use the linearised equations of motion governing small perturbations in the inertial frame. 
These have been discussed and used by many authors (e.g.. Lynden-Bell & Ostriker 1967, 
Papaloizou & Pringle 1978) so we do not present a derivation here writing them in the form: 

j^Z + 2(v • V) + C(|) = VU + i v . (2) 

Here £ is the Lagrangian displacement, C(£) is the linear self-adjoint operator accounting for 
the action of pressure and self-gravity on perturbations (Chandrasekhar, 1964, Lynden-Bell 
& Ostriker 1967). U is the tidal potential given by 

u = gm ( 1 - r ' D(t) ^ ~ 

' - [\r-B(t)\ D{tf )~ D(t) Lm= _ 2>0>2 
Here (r, 9, 0) are spherical coordinates with origin at the center of the planet and polar 
axis coinciding with the rotational axis, all angular momenta being assumed aligned. The 
mass of the star is M, D(t) is the position vector of the star assumed to orbit in the plane 
6 = it/ 2 with D(t) = |D(t)| and $(t) being the associated azimuthal angle at some arbitrary 

* Note that this condition is approximately valid even for a fast rotator like Jupiter. The ratio of Jupiter's rotation frequency 
to Q t is about 0.28 which can be considered as small parameter in perturbation theory (see e.g. Papaloizou & Pringle 1978, 
Lee 1993 for more discussion of this point). 

t An extension of the analysis to take p modes into account is straightforward. 



(3) 
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moment of time t. The moment t — is chosen to correspond to closest approach of star 



and planet. For \m\ = 2, W m = ^J3n/10 and W = —yir/5, with Yj m (0, 0) being the usual 
spherical harmonic. Note, that in equation (3) the indirect term due to the acceleration of 
the coordinate system is included as the second term on the left hand side. On the right 
hand side we take into account only the leading term in the expansion in powers of r/D(t). 

The unperturbed steady state velocity field is v . For the problem considered here, this 
is associated with uniform rotation of the planet such that 

v = ft r x r. (4) 

The viscous force per unit mass is i v . This may be written in component form in Cartesian 
coordinates as 

f: = {i/p)tf = (i/pXMtfl + < - (2/3KO),/?, (5) 

where the comma stands for differentiation and the usual summation convention over the 
Greek indices is adopted from now on. 

The components of the viscous tensor are t a/3 and v is the kinematic viscosity. Since we 
assume uniform rotation of the unperturbed planet, the unperturbed velocity field (4) does 
not contribute to equation (5) so that the velocity components v a may be considered to be 
solely due to the perturbations and are thus related to £ through 

2.2 Normal Modes of the Rotating Planet 

The solution for the perturbations of a planet occurring as a result of an encounter with a 
central star will naturally involve its normal modes £ k with their associated eigenfrequencies 
oo k . These satisfy 

-uUh ~ 2^(v • V)£ fe + C(^) = 0, (7) 

with the standard boundary conditions (e.g. Lynden-Bell & Ostriker 1967). 

When the flow v is absent the eigenf unctions £ fc reduce to £ fc0 and they become those 
appropriate to a spherical star (e.g. Tassoul 1978). The associated eigenvalues are then uj k0 . 
The £ fe0 may then be assumed to be orthogonal in terms of the product (1) and normalised 
by the standard condition 

< Col&o >= !> ( 8 ) 
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In this non rotating limit the eigenfunctions have a simple dependence on <fi through a factor 
exp(im0) with m being the azimuthal mode number. We do not show this dependence 
explicitly here. 

As we use an expansion in the small parameter VL r /VL*\ £ k = £ k o + £,ki and u k = oo ko + uj k i, 
where £ fcl and the frequency correction uj kx are first order in the small parameter. The 
frequency correction has the standard form (e.g. Tassoul 1978; Christensen-Dalsgaard 1998 
and references therein) 

wjfei = -i < &ol( v o- V)& > • (9) 

2.3 Eigenfunction Expansion 

We look for a solution to equation (2) in the form -f 

* = E 6 (*)*€*o(r)- (10) 

k 

The eigenfunctions satisfy the standard equation for a non-rotating object: 

-"lotko + C(* w ) = 0. (11) 

We shall ultimately retain only those corresponding to the I = 2 fundamental mode and 
there are three of these corresponding to azimuthal mode number m = —2, 0, 2. We use the 
standard representation 

Zko = i{r) R Y 2m {9, 0)e r + £(r) s(rVY 2m (9 , 0)). (12) 

The components C(t)r and t;(r)s are real and independent of m, but not independent of 
each other. For a convective isentropic non-rotating planet the motion associated with a 
mode must be circulation free. Accordingly V x £ fc0 = 0, which leads to 

d = kzfr (13) 
dr r 

2.4 Solution of the Tidal Problem 

Substituting equation (10) into equation (2) and using equations (3 - 6) and equations (8 - 
11) with the orthogonality of the eigenfunctions, we obtain under the assumptions of small 
unperturbed rotation rate, small viscous forces and dominance of the / modes: 

b k + u 2 k0 b k + 2iu kl b k + ^(b\ + imVt r b k ) = f k . (14) 

t Adding the corrections £fcl( r ) to £fco( r ) ' n equation (10) does not alter our results. 
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There are three equations of the above form for the modes with to = —2, 0, 2 respectively. 
From now on we accordingly find it convenient to make the equivalence k = to. 
The viscous damping rate 7 has the form 

l = \j ^xvpafal (15) 
where 

<$ = O + & - (2/3)OT- ( 16 ) 
One can readily verify that for a normalised eigenfunction, r j~ 1 has the dimension of time 
and corresponds to a viscous diffusion time across the planet. 
The quantity /J determines the tidal coupling and is given by 

fk =< &>\VU > . (17) 

The coefficients in equation (14) can be obtained in explicit form. Substituting equations (4), 
(11) into equation (14) we obtain (e.g. Christensen-Dalsgaard 1998 and references therein) 

Uki = mf3VL, n (18) 

where 

P=l~ [ Rpl r 2 drp(2^ s + e s )- (19) 
Jo 

An explicit expression for the damping rate 7 is derived in Appendix. The result is 

7 = 4 I** r 2 drpu{^ R - ^ + 3^) 2 + - 3^) 2 + 2& + 3(& - ^) 2 )}, (20) 

Jo o r r r A 

where the prime stands for differentiation with respect to r. Finally, the quantity fl has the 
form 

fl = TW Wme ~ imm > (21) 

where the overlap integral Q is expressed as (PT) 

Q = 2 [ Rpl drpr 3 (U + ^s)- (22) 
Jo 

Equation (14) is close to that obtained by Lai (1997). The only difference relates to the 
treatment of the the last term on the left hand side. We give the explicit expression for 
the damping rate 7 (see equation (20)). Note the contribution proportional to {imVL r b m ) 
in equation (14). This arises from converting Eulerian velocity perturbations to Lagrangian 
displacements. It is necessary to include this term in order to obtain correct expressions for 
the energy and angular momentum exchange between planet and orbit associated with quasi 
static tides (see below). 
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It is convenient to introduce the real quantities 

1 i 
a + = — ^=(£» 2 + 6-2), a_ = -^=(62-6-2), a = 6 , (23) 

and 

/J = ^(tf + /- T 2 ) = V2^W 2 cos (2*(f)), (24) 

/— = -^(fl - £ 2 ) = v/2^^2 sin (2*(f)), (25) 
and rewrite equation (14) in terms of them to obtain 

d + + Ldl a + + aa_ + 7(0+ + 2VL r a_) = /+, (26) 

d_ + cjQ a_ — aa + + 7(d_ — 2f2 r a + ) = /J, (27) 

a + ^oo°o + 7«o = / T ) (28) 
where a = 2\uki I- 



2.5 Energy and Angular Momentum Exchanges 

As we wish to consider tidally induced exchanges of energy and angular momentum between 
planet and orbit, we first consider the energy and angular momentum associated with the 
perturbations. 

The canonical energy E and angular momentum L of perturbations can be expressed as 
(e.g. Friedman and Schutz 1978) 

2' " ' ^ ' ' 2 

and 

a 



E = t{al + al + d 2 ) + ^°(a 2 + a\ + a 2 ), (29) 



L = 2(a + d_ — a_d + — — (a + + a_)). (30) 
Equations governing their time evolution are readily found from equations (26 - 28). Writing 
If = + ^JF' an< ^ ^ = ^ + tne resmt can be expressed in the form 

^ = «o/ T + + ^ /_ T , (31) 
— — = —7(^0 + o? + + a 2 _ — 2£7 r (a + d_ — a_d+)), (32) 

(JjL 

dL T 



dt 
and 
dU 



2{a + f 1 _-a_f 1 + ) (33) 



= -27(a + d_ - a_d + - 2fi r (a 2 + a 2 )). (34) 
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One may regard dE T /dt, dL T /dt, dE u /dt, dh v jdt as determining the increase of the per- 
turbation energy and angular momentum arising from tides and viscosity . These quantities 
can be easily calculated with the help of equations (26 - 28). 

2.6 Quasi-Static Tides 

In the limit when the time-scale associated with changes at closest approach is much longer 
than the internal dynamical time-scale, the planet responds by evolving through a sequence 
of quasi-hydrostatic equilibria. In this approximation, significant perturbations are only 
induced by tides in the vicinity of closest approach. Energy and angular momentum transfer 
to the planet occurs only by the action of viscosity. 

The energy and angular momentum gained by the planet is equal to that dissipated in 
perturbations and can be expressed as 

AE st = - / dt—, (35) 

J — oo CUt 

r+°° dL v 

AL st = - / dt°—, (36) 

J —oo CLt 

where t — corresponds to the moment of closest approach. To evaluate the above integrals, 
we look for the solution to equations (26 - 28) in the form of an expansion in terms of the 
small parameter f^/f^, where 



GM 



mm 



and D m i n is the value of periastron distance. 

We also assume that u 00 ~ and Q r ~ Q p . In the leading approximation the solution 
is trivial: 

f T f T 

ao = — j-, a± = — g-. (38) 
w oo ^oo 

Substituting (38) into (32) and (34), and evaluating the integrals in (35) and (36) we obtain 
AE st = -^-(h - 2Q r I 2 ), AL st = ^(I 2 - 20J 3 ), (39) 

^00 ^00 

where 

h = r dm?) 2 + uif + (/. T ) 2 ), h = r dt{fj T _ - an 

J —oo J —oo 

h= r°dt(( A? +(£?). (40 

J —oo 

For a parabolic orbit, the evaluation of the integrals in equations (40) is straightforward. 
After elementary, but rather tedious calculations we obtain 
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AEx=Pl 7^^9l (1 _m^ (41) 

^oo D min 117v2n p 
AL^A^^d-l-^-), (42) 

^00 ^min 00 V ^ 'p 

where fi 1 = ^^■y/2n 2 w 31,6, /3 2 = §§§7r 2 21,4. It is very helpful to introduce new 
dimensionless variables in equations (41 - 42). Instead of using D min directly, we specify the 
closest approach distance by use of the Press- Teukolsky parameter 77 defined as 



v = 77- 



m pi D m i n (43) 



n p \m i$ ' 

The linear theory of perturbations is valid only for sufficiently large values of 77. As we will 
see in the next Section, in our problem we have typically 77 ~ 10, and the use of the linear 
theory is fully justified. 

We find it convenient to express Q, u 00 , 7, Q r in terms of natural units. In so doing 
we introduce dimensionless variables indicated with a tilde. These are specified through 
Q = Q\fm pl R p i, uj 00 = u 00 Q*, 7 = 7^*, tt r = fi r ^*- Note that 7 _1 , being related to 
the internal viscous diffusion time is expressed in units of the internal dynamical time. 
Thus although it is dimensionless, 7 may be a very small quantity. Using the dimensionless 
variables, the energy and angular momentum changes may be expressed as 

AF tf^^M U2 ^M f44^ 

AEst ~ (44) 

and 

AL st =^^(l-A. (45) 
Here these are normalised in terms of the energy and angular momentum scales E p i = 



GrripjRpi and L p i = m p i ^Gm p iR p i respectively. However, these values may only be ap- 
proached for very close encounters and comparable viscous and dynamical times. 

Equations (44) and (45) are formally identical to the equations obtained by Hut (1981) 
who used a weak friction model in the limit of the planet orbit being marginally bound to 
the star, provided that we make the identification 

i = (46) 

where k is the apsidal motion constant, T = l/(Jl*r), with r being the time lag between the 
tidal forcing and the response of the planet that occurs through frictional processes. 

The advantage of our approach is that it allows us to obtain the correct expressions for 
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the energy and angular momentum exchange due to the stationary tides. In this way the 
otherwise unspecified time lag r may be expressed in terms of quantities determined by the 
form of viscosity in the planet, the planet's structure, and also by the mode of oscillation 
responding to the tidal field. Note that in framework of the standard theory of quasi-static 
tides the similar relations have been obtained by Zahn (Zahn 1989 and references therein) 
and by Eggleton, Kiseleva and Hut (1998). 

As can be seen from equations (44) and (45), in the limiting case rj — > oo, the planetary 
angular velocity tends to evolve faster than its orbital energy with the consequence that if 
there are repeated encounters, the planet achieves a state of pseudo-synchronisation with 
f2 r = il ps = <C fi*, at which point, exchange of angular momentum between the 

orbit and the planet ceases (Hut 1981). For £l r = n ps the energy gain follows from equation 



where (3^ = jg^Pi ~ 6.64. We comment that the angular velocity f2 ps roughly corresponds 
to the orbital angular velocity at closest approach. In the next Section we use equation (47) 
to estimate the importance of stationary tides for the orbital evolution of the planet. 

2.7 Dynamic Tides and the Tidal Excitation of Normal Modes 

When the tidal encounter with the central star occurs on a time-scale which is long but 
not extremely long compared to the internal dynamical time the planet can be left with a 
normal mode excited at some non-negligible amplitude. There will be energy and angular 
momentum associated with this mode that is transferred from the orbit. In contrast to the 
situation for quasi-static tides the effect does not depend on viscosity. We refer to it as the 
non- stationary or dynamic contribution to the energy and angular momentum transfer. 

In order to calculate the non-stationary contribution to the energy and angular momen- 
tum gain we can neglect the viscosity term oc 7 in equations (26 - 28) which may then be 
written in the form 

a± + ujQ a± ± aa T = /J, d + u% a = fo- (48) 

We use below the solution to equation (48) obtained by the method of variation of parameters 
correct to first order in cr/c^oo in the form 

\ ft / , , , , 

a± = / dt {f±(smuj±(t — t ) + sin uj t (t — t )) + f T (cosuj±(t — t ) — cos LU T (t — t ))},(49) 



(44) as 



AE st = fo^ 




(47) 



2cj o 



—00 
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a = / dtf smuj 00 (t - t'), (50) 

where cu± = c^oo±f • In equations (49 - 50) we assume that the perturbations of the planet are 
absent before the close encounter (t = — oo) (the so-called first passage problem). General 
expressions appropriate to repeated encounters may easily obtained from the expressions 
derived for the first passage problem (see below). 

In the limit t — > oo the coefficients multiplying sin(u;oo(i)), sin(u;±(i)), and cos(u;oo(t)), 
cos(u±(t)) in equations (49 - 50) contain converging integrals. These contain the informa- 
tion about the residual mode excitation after the encounter. Extending the upper limit of 
integration in these integrals to infinity, we have 

a + = — !— sin(cj + t) + I-(u)-) sm(u-t)), (51) 

2^oo 

a_ = (J_(u;_) cos(u^t) — I + (oj + ) cos(uj + t)), (52) 

2^oo 

and 

a = — I (oj 00 ) sm(u 00 t), (53) 
^oo 

where 

/+oo f+OO 

dt(f + cos(c<jt) ± /_ sin(atf)), I (uj) = / dtf cos(ojt). (54) 

-oo J— oo 

These integrals can be represented in the form 

/ ± = 2V2,^-K m , I = 2,^-K , (55) 
where m = =f2 for (±), and 

K m = ^ f°° cosM + m<t>(0), (56) 

Z7T J-oo X[t) 6 

and x(t) = D(t)/D min . The integrals (56) have been calculated numerically and approxi- 
mated analytically by Press and Teukolsky (1977). Lai (1997) has obtained useful analytical 
expressions for these integrals in the limit (our]) — > oo in the form 



x-^^jvr^-^ K ^-^7W {l+ tr> K *«fti-w*-m^ 

where z = y/2urj. In that limit we have |if_ 2 (2)| > |-Ko(z)| > 1-^2(^)1 

Substituting equations (51 - 52) into equation (30), we obtain the gain of angular mo- 
mentum after "the first" periastron passage 

AL ns = J-{/> + )-/ 2 (c,_)}. (58) 
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Analogously substituting equations (51 - 53) into equation (29), we obtain the expression 
for the gain of energy 

where in the last equality we distinguish between u>± and c^oo only in the factors depending 
on uj exponentially. 

In contrast to related expressions calculated in the approximation of quasi-static tides, 
the expressions (58 - 59) give the energy and angular momentum gain associated with some 
particular (in our case, fundamental) mode of pulsation. The transfer of the energy and 
angular momentum from the mode to the planet proceeds during some dissipation time 
tdiss = Too 1 which could be much larger than the orbital period of the planet. Note that, 
in general, the mode damping rate 700 is different from the damping rate 7 introduced for 
quasi-static tides. 

The expressions (58 - 59) can be significantly simplified. For a non-rotating planet (Q r = 
0) and large values of rj, we can neglect the contribution of (— ) and (0) terms, and obtain 
from equations (55), (57) and (59) 

AE ns (Q r = 0) - - ^^Q^e-^^E^ (60) 

and AL ns m ~ 2 ^" s ■ Of course, these expressions follow directly from the general 

expression for the energy gain provided by PT and the general relation between the mode 
energy and angular momentum for m = — 2 (e.g. Friedman and Schutz 1978). 

In the general case of a rotating planet, we can substitute the asymptotic expressions 
(57) in equation (58) and obtain 

AL ns « ^(^ooQ) V^ } (1 " je^)L pl , (61) 

15 2 14 (^oo^) 

where a = a/VL*. 

It follows from this expression that the non-stationary tides can spin up the planet up 
to the angular velocity which when inserted into the right hand side of equation (61) gives 
no additional increase in planet angular momentum or AL ns = 0. It is readily seen that this 
angular velocity is given by 

= = 7-7=^ ln (il^ooV) « ^ In (6.53^), (62) 
where we have used equation (18) and the definition of a. 
Typically (5 ~ 0.5, uqq ~ 1 and we have 
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2 + In r\ 

^crit ~ (63) 

V 

It is interesting to note that for distant encounters CoqqT] may be much larger than unity. This 
has the effect that, in contrast to the situation with quasi-static tides, the critical angular 
velocity may easily exceed the orbital angular velocity at periastron. For closer encounters 
with sufficiently small values of r\ (but still larger than one) the critical angular velocity 
can even be of the order of the characteristic frequency Q*. This comes about because in 
simple terms, the pattern speed of the excited normal mode is also characteristically and 
the star has to rotate by an amount approaching this before the mode excitation process is 
significantly affected. 

Thus in principle, the non-stationary tides might lead to rotational break-up of the 
planet at values of periastron beyond the Roche limit. However, our approximations are 
clearly invalid at such high angular velocities. 

Substituting equation (57) into equation (59), we have the similar expression for the 
energy gain 

AE n , » i^QV-^^d + — + —°e*<»>)B+ (64) 
15 2 6 (c^oo^) 2 14 (cj o^) 

It is easy to see that the energy gain attain its minimal value at Q = Vt crit . Taking into 
account that the expression in the round brackets in equation (64) is equal to 4 for Q = Q C rit, 
we have 

AE mm = AE ns (Q = n cnt ) = 1-^% = -L^-e-^^E^. (65) 

32 {Ufi]) 5y2 V 

2.8 Energy and angular momentum gain resulting from multiple encounters 

A planet in a highly eccentric orbit will undergo multiple periastron passages each of which 
is associated with an energy and angular momentum exchange. If dissipative processes are 
sufficiently weak, normal modes remain excited between periastron passages with the conse- 
quence that successive passages are not independent. We here analyse such a situation and 
find that for highly elongated orbits there may be a stochastic build up of normal mode 
energy at the expense of orbital energy and provide an estimate of when this will occur. 

As pointed out by several authors (PT, Kochanek 1992, Mardling 1995 a,b) the amount 
of energy and angular momentum gained (or released) during a close approach by some 
particular mode of oscillation depends on the mode amplitude and the phase of oscillation 
at the moment of periastron passage. Let us consider this effect in more detail. As indicated 
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above, we assume that the mode damping time ~ 7J 1 is much larger than the orbital period 
of the planet 



Porb = l^L (66) 

and thus neglect the effect of dissipation. Here E orb = is the specific binding energy of 
the planet's orbit (from now on called the orbital energy) and a or b is the semi-major axis. 

We assume that the tidal interaction occurs impulsively. The energy and angular mo- 
mentum exchanges are assumed to occur at the moment of periastron passage and at any 
other time the mode evolution is described by solutions of equations (26-28) without external 
forcing terms in the form 

a + = B + {j) sin(^ + t + Mj)) + B.{j) sin(w-tO') + 1>-U)), ( 67 ) 
a_ = B_(j) cos(cu_t + V-0')) - B + {j) cos(cM + ^ + (j)), (68) 
a = V2B (j) sin(cu 00 t + Mj))- (69) 

Here Bi(j), ipi(j) are the mode amplitudes and phases after (j — 1) periastron passages and 
we have introduced the index (i) which can be one of (+, — , 0). 

The phases have been adjusted to make the time t — correspond to the jth periastron 
passage. To take into account the interaction at the jth periastron passage we appropriately 
add in the expressions (51-53) to (67-69). 

The result can also be represented in the generic form (67-69) but with new amplitudes 
Bi(j + 1) and phases ipi(j + 1) given by 



B,(j + 1) = JW)+% + 2B,m<M*iW), (7°) 

^ + »-^ mSt))U )+m ' (71) 

where 6< = ^ for % = (+, -) and b = (see also Mardling 1995 a, Mardling & 

Aarseth 2001). The additional phase <f>i(j) = uJiP or b(j) being the normal mode phase shift 
during the orbital period immediately after the jth periastron passage is added in to make 
the time t — now correspond to the j + 1th periastron passage. We assign the index (j) 
to the orbital period stressing the fact that the period can change with time. Obviously, the 
first term in (71) ~ when B^j) 3> foj. 

The expressions for the mode energy and angular momentum just before the jth perias- 
tron passage follow from equations (29-30) together with equations (67-69) as 
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Em(j) = E ^t^ BfU) - -o 2 o E BfV), (72) 
L m (j) = 2u 0Q (BP + V)-BlV)). (73) 

The gain of energy and angular momentum after the jth periastron passage follows from 
equations (70), (72) and (73) as 

AE m (j) = u& E ((&? + 2 W) cos(^O')))) , (74) 

i 

AL m (j) = 2LO 00 {b 2 + + 2b + B + (j) cos(0 + (j)) - (b 2 _ + 26_ J B_(j) cos(0_(j)))}. (75) 

According to equations (74) and (75), the energy and angular momentum can either increase 
or decrease after periastron passage depending on sign of cos(ipi(j)) and the value of Bi(j) 
(PT; Kochanek 1992; Mardling 1995 a,b). 

2.9 A two dimensional map for successive encounters 

Equations (70) and (71) define an iterative map. This can be brought into a simpler form 
by introducing the two-dimensional vectors Xj(j) with components 

x\ = cos(Mj)), x? = ^ MM))- (76) 
bi bi 

We obtain from equations (70) and (71) 

x i O' + l)=R(0 i )(x i O')+e), (77) 

where R(0,) is a two-dimensional rotational matrix corresponding to a rotation through an 
angle —<f>i, and e is a two-dimensional constant vector with components e 1 = 1, e 2 = 0. The 
mode energy is expressed in terms of the new variables as 

£ m (j) = A£„ s EA|x,| 2 , (78) 

i 

where AE ns is given by equation (64), and the dimensionless coefficients $ = ^§r^-. Note 
that EiA = 1- 

For a constant period P or b{j) = -Porb, the angles 0j are constant, and the general solution 
for the mapping given by equation (77) can be easily found in the form 

XiO') = xf + A(j'0,)x o , (79) 

where the constant vectors xf give the stationary point under the mapping defined by 

equation (77) and are found from the condition that x;(j + 1) = Xj(j). They have the 
components 
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st(l) 



1 

2' 



st(2) 



sin(0j) 



(80) 



2(1 - cos(&))" 



A(j0j) is the two-dimensional rotational matrix corresponding to a rotation through an 
angle —j4>i, and the constant vectors x° are determined by the initial conditions. According 
to equation (79), the vectors Xj(j) lie on closed circles centred on the stationary points x s * 
with radii determined by the initial conditions. Obviously, no persistent transfer of energy 
from the orbit to the planet is possible in that case. 

When the orbital period changes, the situation is much more complicated. It is clear from 
equations (76) and (77) that only the values of (fii(j), mod 2tt are significant. Provided 
that the change of orbital period from one periastron passage to the next is sufficiently large 
(namely AP or b > 2-k/ujw), successive phases 0«(j) mod 2tt, and consequently the angles 
ipi(j) behave in a complicated manner. These quantities can be approximately represented as 
uniformly distributed random variables. With this assumption we can average over cos(?/>i(j)) 
in equations (70) and (74) to obtain 



This indicates a secular increase of mode energy on average and thus may be thought 
of as a stochastic instability. Note that contrary to the usual random walk problem, in 



our case dispersion grows as fast as the average value: y< E^ > — < E m > 2 ^< E m >. 
Nevertheless, in our estimates of orbital evolution of the planet we assume that the energy 
gain is determined by from equations (81). 

In a standard approach to the problem the change of the period is caused by exchange 
of energy between the orbit and the mode due to tidal interaction. Assuming that there is 
no dissipation and radiation of energy away from the planet, the sum of the orbital energy 
and the mode energy is conserved, and the Keplerian relation (66) can be used to express 
the change of orbital energy in terms of change of the mode energy and close the set of 
equations (77) and (78). Mardling (1995a, 1995b) numerically investigated the problem of 
evolution of the orbital parameters and the mode energy for two convective polytropic stars 
of comparable masses. It was found that systems with a substantial eccentricity exhibit 
stochastic behaviour. The origin of this result is clear from qualitative point of view. Indeed, 
more eccentric orbits have larger semi-major axes a, and therefore larger fractional changes 
of the orbital period following from a periastron passage. We have 



< B t {jf >= b 2 J, < E m >= AE ns3 . 



(81) 




AP orb (j) 



3 P orb AEjJ) 

2 Eorb ™>pl 



= —6n 




a 5 AEJJ) 
mf m pl 



(82) 
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For our problem it is important to establish a condition of onset of stochastic instability. 
To do this we consider a situation where the mode energy is much smaller than the orbital 
energy of the planet, so that the orbital period does not change significantly from some 
'initial' value P . In that case we can express the angles fa in terms of the mode energy 
E m (j) through 

fa(j) = fa(0)-a l E m (j), (83) 

where fa(0) = UiP . The dimensionless quantities determine how fast fa(j) changes with 
E m (j), and follow from equation (82). Thus 
Uja 5 / 2 AE ns 

Cti = 6ll rrr . (84) 

{GM f 2 m pl 
The dimensionless mode energy E m (j) is defined as 

E m {j) = E m (j)/AE ns = £AN 2 . (85) 

i 

Equations (77) and (83 - 85) form a closed set of equations, which depend parametrically 
on the values of fa(0) mod 2ti, c*i and As initial conditions we can use the result of the 
first passage: x}(0) = 1, xf(0) = 0. 

We analyse the evolution of this set of equations for Q r = 0. In that case we can suppress 
the index % in all expressions and the problem depends parametrically only on the values 
of a and 0o mod 2n. For a fixed value of 0o and sufficiently small values of a, the solution 
generated from equation (77) is attracted to a stationary point near to the one that exists 
for a = and which satisfies equation (80) (see figure 1). We comment that the map is not 
area preserving in general and that for a ^ 0, there are an infinite number of fixed points 
that have angles that can be found from the equation 

* = * (0) - 2(1-1(0)) - (86) 
When a exceeds some value a*, this stationary point described above becomes unstable, 

and the solution exhibits a complicated behaviour, but stochastic growth is not observed 
(see figure 2). Finally, when a exceeds some critical value a cr u, the solution starts to exhibit 
stochastic behaviour and the mode energy grows on average (see figure 3). Analytical cal- 
culation of a crit is a very complicated mathematical problem so we use a simple numerical 
method in order to establish the dependence of a crit on fa. Namely, we iterate the map 
(77) for a large number n of iterations (n — 5 • 10 4 , 10 5 , 1.5 • 10 5 ) and different values of a 
and fa and find the boundary between stochastic and stable solutions. We solve the system 
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Figure 1. The evolution of the mapping (77) is shown in the (a; 1 , a; 2 ) plane for a = and a = 0.1. Circles correspond to 
positions of the vector x(J) for different j. Circles connected by the solid and dotted lines represent the cases a = 0, a = 0.1, 
respectively. The angle <po = n/2. Both curves start at the initial value of x(0) with coordinates: xj(0) = 1, £ 2 (0) = 0. The 
closed curve corresponding to a = is a parallelogram where successive values of x(j) are situated at vertexes. The curve 
corresponding to a = 0.1 is a spiral which ultimately approaches a stationary point with coordinates defined by equation (80). 




Figure 2. Same as figure 1, but with a = 0.3. We represent here 5000 point corresponding to 5000 iterations of the map (77). 
The stationary point referred to in the text is unstable in this case, but all points are confined to a small region around it. We 
have checked for a much larger number of iterations that |x(j)| remains bounded. The case is intermediate between attraction 
to a fixed point and stochastic behaviour. 



starting with a = 0, and steadily increase a until E m {n) exceeds the iteration number n. 
The corresponding value of a is associated with a crit . We present the dependence of a crit 
on 0o i n figure 4. One can see from this plot that this value very weakly depends on the 
iteration number n. For our estimates below we use the averaged value of a crit 
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200 



Figure 3. Same as Fig. 1 and Fig. 2, but with a = 0.35. This value of a approximately corresponds to a cr u for 0o = tt/2 (see 
Fig. 4). The evolution can be approximately described as a jumping of x(j) for successive j between different points on circles. 
The radius of the circles (proportional to the mode energy) grows on average. 



3 i 




Figure 4. The dependence of a cr u on </>q is shown. The solid, dashed, and dotted lines correspond to n = 1.5- 10 5 , 10 5 , 5- 10 4 , 
respectively. 



1 f 2n 

dcrit = — / d(p a crit w 0.67. (87) 
2n Jo 

One can see from equation (84) that a oc a 5 / 2 . Therefore, the stochastic instability always 
sets in when the semi-major axis a is sufficiently large: a > a st = a(a = a crit ). To calculate 
a st explicitly one should specify a particular model of the planet. We calculate it in Section 
3.5 below. 
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The case of rotating planet can be treated by a similar way, but in this case we have a 
large number of parameters defining the map. The increased number of degrees of freedom 
in that case is likely to result in decrease of a cr u with Q r . Therefore, we assume that we can 
use a cr it(fl r = 0) as an approximate upper boundary for the rotating case. 

3 ORBITAL EVOLUTION OF A MASSIVE FULLY CONVECTIVE 
PLANET DUE TO TIDAL INTERACTION 

In this Section we apply the results of the previous section to estimate the effect of tides 
induced by the central star on the orbital evolution of a massive planet. Consider a planetary 
system consisting of several massive fully convective planets orbiting a central star at typical 
separation distances ~ 1 — lOOAu. Suppose that one of these has lost a significant amount of 
angular momentum (most probably due to interaction with other planets) and has settled 
into a highly eccentric orbit. The orbital semi- major axis a would then be of the order of 
'initial' separation distance ~ 1 — lOOAu. The minimal separation distance D min is related 

L 2 

to the orbital angular momentum per unit of mass L orb through D min = and could 
be of the order of several stellar radii R*. Situations of this type occur in simulations of 
gravitating interacting planets (Papaloizou & Terquem 2001, Adams & Laughlin 2003). If 
close enough approaches to the central star occur, tidal interaction could become important 
for determining the subsequent orbital evolution of the planet. Below we estimate the relative 
contribution of stationary and dynamical tides, both exerted by star on planet and planet 
on star on the orbital evolution. In particular we specify at which separation distances we 
can expect tidal effects to become significant. 

If we assume that the orbital evolution of the planet is determined by tides, while the 
orbit remains highly eccentric, the total relative change of the orbital angular momentum is 
much smaller than that of the orbital energy. Thus we may suppose the orbital momentum to 
be approximately conserved. That leads to the well known result that the 'final' separation 
distance of the planet in the 'final' circular orbit attained at the end of tidal evolution, Df in , 
is twice as large as D min . The shortest observed periods are typically 3 days for Jovian mass 
planets around solar type stars (Butler et al 2002). Thus D m i n is directly related to the 
observed period P Q b s of the planet through 




(88) 
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where R Q = 7 ■ 10 10 cm is the radius of the Sun, and P3 = P bs/ (Sdays). Using equation (43) 
we can also express the parameter r\ in terms of P3: 



For a planet of one Jupiter mass and radius, and a central solar type star, 770 ~ 1- 
3.1 The planet model 

For a close encounter with the central star occurring with a given distance of closest ap- 
proach, the strength of tidal effects acting in a planet is sensitive to the magnitude of its 
radius. This is because the energy gain through the action of dynamical tides decreases ex- 
ponentially with i] (see equation(64)) and rj oc R pl 3 ^ 2 is a monotonically decreasing function 
of the radius of the planet. The energy gain associated with stationary tides being oc i]~ 5 
also decreases rapidly with 77. 

In order to determine the radius, a model for the structure of the planet taking account 
its history is needed. The radius decreases as the planet cools on a thermal time-scale. 
Accordingly the model should specify the planetary structure and radius as a function of 
age. 

A young gravitationally contracting planet will be fully convective. If we assume that 
a turbulent viscosity arising from convection determines the damping rate 7, it will be a 
function of the planet radius and luminosity both of which are functions of time. In this 
case too, estimation of the importance of tidal interaction requires specification of a time 
dependent model of the planet. 

Because convection is extremely efficient, a fully convective planet can be considered 
as an isentropic body. § For a given composition isentropic models of a given mass can 
be conveniently parametrised by the planet radius R p i. We have calculated such models 
using the equation of state of Saumon, Chabrier & Van Horn (1995) for solar composition 
neglecting the possible presence of a central small solid core. We have considered the masses 
m p i = 1Mj,5Mj, where Mj ps 1.9 • 10 30 g is the mass of Jupiter, with radii in the range 

§ Note that this approximation may break down at a boundary between the molecular and metallic phases of Hydrogen and 
Helium if the phase transition is of the first kind. See also Discussion. 




(89) 



where 



m p i Rl 



(90) 
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Figure 5. The dimensionless density p as a function of the dimensionless radius r/R p i. The three solid lines represent m p i = 
lMj with R p i = IRj, 1.5-Rj and 2Rj respectively. The three dotted lines are for m p i = 5Afj with R p i = 1.03i?j, 1.54i?j and 
2.06-Rj. The central density grows with R p i for both masses. 

IRj — 2.1Rj, where the Jovian radius Rj ~ 7- 10 9 cm. The dimensionless density distribution 
p(r) as a function of the dimensionless radius r/R p i is shown in figure 5. Note that the models 
with larger radii are more centrally condensed, the effect being less prominent for m p i = 5Mj. 
Also note the sharp change of the density at r/R p i « 0.82 for the model with m p i = lMj and 
R p i = IRj. This is due to the phase transition between molecular and metallic Hydrogen 
that has been assumed to be first order. This effect is absent for the models with larger radii 
and larger mass. For these models, the corresponding adiabats have larger entropy and go 
above the critical point of the Saumon, Chabrier & Van Horn (1995) equation of state. We 
calculate the dimensionless frequencies of the fundamental mode u) o and the dimensionless 
overlap integrals Q for these models. The results of these calculations are presented in figures 
6 and 7. 

The dimensionless frequencies increase, and the dimensionless overlap integrals decrease 
with R p i these effects being less prominent for m p i = 5Mj. The value of uj 00 for the model 
with mpi = lMj, Rpi = IRj is in a good agreement with earlier calculations by Vorontsov 
(1984), Vorontsov et al (1989), Lee (1993). 

The coefficients (3 (equation (19)) determining the frequency splitting due to rotation 
are close to 0.5 for all models. Note that this would be exact for an incompressible model. 
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Figure 6. The dimcnsionless frequency ujoo as a function of the radius of the planet (in units of the Jovian radius). The solid 
and dashed curves correspond to m p i = lMj and 5Mj respectively. 




Figure 7. The dimcnsionless overlap integral Q as a function of the radius of the planet (in units of the Jovian radius). The 
solid and dashed curves correspond to m pi = lMj and 5Mj, respectively. 

3.1.1 Evolution of the radius with time 

In order to find the dependence of the planetary radius on time, a model of the planetary 
atmosphere is needed. For a given mass and radius, this is required to determine the lumi- 
nosity which in turn determines the cooling rate as a function of radius. This is a highly 
non trivial problem given significant difficulties in calculating molecular opacities for the low 
temperatures and high densities of interest. Other complications may arise from the possi- 
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Figure 8. The radius of the planet (in units of the Jovian radius) as a function of time (in units 10 6 j/r). The solid and dashed 
curves correspond to m p ; = lMj, 5Mj, respectively. Note that the dependence of radius on time is rather similar for both 
models. 



le+31 



le+30 




le+24 



Figure 9. The internal luminosity of the planet (erg/s) as a function of time (in units 10 6 j/r). The solid and dashed curves 
correspond to m p i = lMj and 5Mj, respectively. The more massive planet is more than ten times more luminous. 

ble presence of clouds of different condensed chemical species, etc. We use the relationship 
between radius and luminosity and hence the planetary age calculated by Burrows et al 
(1997) in the framework of a non-gray atmospheric model. The results of these calculations 
are shown in figures 8 and 9. As one can see from figure 8, the radius of the planet decreases 
from ~ l.TRj at t = 10 6 yr to ~ Rj at t ~ 8 • 10 9 yr. Note that this result is close to that 
obtained by Graboske et al (1975) for their model of the evolution of Jupiter. We here point 
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Figure 10. The convective time scale t zh as a function of time (in units 10 6 yr). The solid and dashed curves correspond to 
m p i = lMj and 5Mj respectively. Note that this time scale is the same order of magnitude for both masses. 

out that tidal heating might significantly delay or even reverse the cooling of the planet. 
This effect is not considered here as we focus on the initial stages of tidal evolution and so 
do not solve a self-consistent problem for the thermodynamics of a planet with an energy 
source due to tidal heating. 



3.1.2 The convective time scale 

An important time scale is the characteristic time associated with convective eddies t c . These 
will be of large scale in the fully convective planetary interior. This characteristic time has 
been estimated by Zahn (1977) to be 

t c ~^Pf^) 1/3 , (91) 

where E th is the energy transport rate due to convection. This is comparable to the total 
luminosity. The dependence of t z h on time is shown in figure 10. 

One can see from figure 10 that t z h grows with time from ~ 0.5yr at an age of 10 6 yr to 
~ lOyr at an age of 10 9 'yr. Assuming that dissipation is determined by a turbulent viscosity 
that acts like an anomalously large molecular viscosity, we can express the dissipation rate 
given by equation (20) as 

7 = X ( m pii R pi) ? ( 92 ) 

tzh 

where the coefficient X(m p i, R p i) ~ 10~ 2 — lO^ 1 for different models of the planet. 
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In fact the dissipation rate may be much smaller than that given by equation (92). This is 
because as the largest convective eddies have the size of order of the mixing length, which in 
the deep interior is on the order of the local radius, they can only fully distribute the energy 
and momentum of the hydrodynamical motion induced by the tidal field on the mixing 
length scale in a time ~ t c ^> tf, where tf ~ lday is the characteristic time associated 
with tidal forcing (Goldreich & Nicholson 1977). This leads to a correction factor ~ {^-) n in 
equation (92). In the literature it has been suggested that n may be equal either to 2 (e.g. 
Goldreich & Nicholson 1977) or to 1 (Zahn 1989). Neglecting the contribution of layers close 
to the surface of the planet, we adopt t c ~ t z h, and modify the expression (92) taking into 
account the correction factor so obtaining 

i m - ( ^H A(m r Rpi) ). 03) 

tzh *zh 

3.2 Orbital evolution due to quasi-static tides 

We show here that in our approximation, the time scale for orbital evolution due to quasi- 
static tides of a planet in a highly eccentric orbit is exceedingly large. In order to make an 
estimate of the evolution time scale we assume the planet to be rotating at its 'equilibrium' 
angular velocity Q ps and therefore evolving at fixed orbital angular momentum. This is 
reasonable on account of the low moment of inertia associated with the planetary rotation 
compared to that of the orbit. The planet in a highly eccentric orbit can be considered to 
have a sequence of impulsive energy gains occurring at periastron passage. Each energy gain 
AE st is given by equation (47). Since the energy gain depends on orbital parameters only 
through the dependence of 77 on the minimum separation distance and for highly eccentric 
orbits the minimum separation distance is determined by the value of the orbital angular 
momentum, the energy gain AE st is constant during the orbital evolution. 

The law of energy conservation gives 
AE st GMm„i . . 

ivf + ^ o = ' (94) 

where a is the semi-major axis, and P or b = 2ir^ is the orbital period (also see equation 
(66)). From equation (94) we obtain the evolution equation for the semi-major axis in the 
form 

- = - — , (95) 

(X b eV 

with 
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(GMm pl \ I a 3 



It is convenient to introduce a scale a* defined by the condition that the planetary internal 
energy is approximately equal to the orbital energy in the form 

^i^M.NfflW.. (97) 
m pl \RjJ \M Q m pl J 

Using this we can express t ev in the form 

E p i /a' 
2AE st W 



where 



M \ ( MjR 



One can easily estimate from equation (98) that t ev is very large. Indeed, AE st /E p i can 
be estimated from equation (47) to be 

AE st /E pl ~ 1 (100) 

where we have assumed that cj o and Q ~ 1. Taking into account that ~ 1 — lOyr, and 
A ~ 10~ 2 , we find 7 using equation (92) to be 7 < 10~ 6 . Since rj ~ 10 in our problem, we 
obtain AE st /E p i < 10 -11 , and thus 

a 1/2. Mj, 3/2 
~ 10 11 - ^ yr. 
a* m p i 

More accurately, we can take into account the fact that AE st depends on several quantities 
which are functions of time, making in turn t ev a function of time also. We show the depen- 
dence of t ev on t in figure 11 assuming that f] = 10, a = a* and the damping rate is given 
by equation (92). One can see from figure 11 that the time scale of orbital evolution due to 
quasi-static tides is very large, being typically the order of ~ 10 9 — 10 10 yr for m pl = 5Mj, and 
order of magnitude larger for m p i = Mj. In fact, the evolution time may even be much larger 
if we take into account the correction factor to the damping rate (equation (93)). Therefore, 
if we assume that the viscous damping rate is determined by convective turbulent viscos- 
ity, the evolution of a planet's orbit due to quasi-static tides appears to be impossible for 
semi-major axes <: 0.1 Au. 

We have also made simple estimates of the evolution due to quasi-static tides acting in 
a star with a solar-like convective envelope and found that the rate of evolution of a planet 
orbit with orbital parameters typical of our problem due to tides in the star also seems to 
be rather weak. 



3/2 
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10000 



Figure 11. The time scale of evolution due to quasi-static tides, t ev , (in units of 10 6 yr) as a function of time ( in the same 
units) for a = a*. The solid and dashed curves correspond to m p i = lMj and 5Mj respectively. 

3.3 Dynamic Tides in the star versus Dynamic Tides in the planet 

In order to estimate effects due to dynamic tides in the star we model it as a simple n = 3 
polytrope and calculate the energy gain due to dynamic tides, AE ns * using the expressions 
given by PT. Our results agree well with those of Lee and Ostriker (1986), but we include the 
contribution of a larger number of g modes to AE ns * and therefore obtain a slightly larger 
value of AE ns * for large separation distances. It is important to note that the quantity i]* 
associated with tides in the star is different from that associated with the planet (equation 
(43)) being given by the expression (PT) 



rj* 



n 3/2 d 3/2 

- J-' 777.7,77. \ ^ / ^0 \ 



9(nr) P* 
-tt* 



;ioi) 



In order to show the values of AE ns and AE ns * in the same figure, we use the orbital period 
P b s as independent variable. These quantities are represented in figure 12 and figure 13 for 
rripi = lMj and m p i = 5Mj respectively. One can see from figure 12 that the energy gain 
AE ns sharply decreases with P b s - However, in most cases AE ns 3> AE ns * at P t, s ~ lday. For 
a non-rotating planet with m pl = lMj tides in the planet are dominant for all interesting 
values of AE ns : AE ns /E p i > 1CT 10 — 1CT 9 . Even when the planet rotates at = Vt crit 
(corresponding to minimal possible energy gain for a given P b s , see equation (65)), the tides 
in the planet dominate for models with R pi > lARj. The planet with m p i = 5Mj is more 
dense, and the role of tides operating in the planet is less prominent. In that case the curves 
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P obs (days) 



Figure 12. The energy gains AE ns , AE„ st obtained by the planet (solid and dotted curves) and by the star (the dashed curve) 
after 'the first' periastron fly-by. The six solid curves correspond to a non-rotating planet (Q r = 0) and represent the expression 
59) for the six different values R p i = lRj, 1.2Rj, lARj, l.GRj , l.8Rj and 2Rj. Obviously, the curves corresponding to larger 
values of Rj have larger values of AE ns . The six dotted curves represent the same quantities but calculated for a planet with 
'critical' rotation rate Q r = Q C rit- They are given by the expression (65). 




Figure 13. Same as 12, but calculated for m p i = 5Mj. 

corresponding to 'critical' rotation are situated below the curve corresponding to the tides 
in the star for almost the whole range of P Q b s - 

3.4 Orbital evolution due to Dynamic Tides 

If we assume that either each impulsive energy transfer to the planet is dissipated before the 
next, or that the condition of the stochastic build-up of mode energy is fulfilled, in either 
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case we can write the equation for the evolution of the semi-major axis due to dynamical 
tides in a form similar to equation (95) 

hl ° 1 (102) 



a w tio(t)y/aio ' 
where aio = a/10 Au. The time 

is the 'local' evolution time for the orbit with semi-major axis ~ lOAu ^> a* and e ns = 
10 9 (AE ns + AE nsif )/E p i. We would like to stress that if the tidal energy transfer process is 
in the stochastic regime, equation (102) must be understood as describing the evolution of 
an average value of a (see previous Section). 
Equation(102) can be integrated to give 



dt' 



aio = a in I 1 - —— I ——7- ] = a in I 1 - ] , (104) 

\ 2^fa~ J t in t 10 {t ) J \ \ a in J 

where aj n and tj„ are the 'initial' values of the semi-major axis (in units of lOAu) and the 
initial time respectively. The quantity 

a ^\(L^) < 105 » 

defines the typical 'capture' scale. This is because at any moment of time t, the orbits 
with initial values of the semi-major axes ai n < a cap (t) have changed their semi-major axes 
significantly before the time t and the orbits with a in > a cap (t) have a ~ a in . 

We can define the evolution time more accurately as the time when the term in brackets 
in equation (104) is equal to zero. Thus 
1 /•*«> dt' 

V7^ / T17\ = L ( 106 ) 
Obviously, this time is a function of ai n and t{ n . For a constant value of tio, t ev ps 2 v /a^ti . 

The results of calculation of t ev and a cap are presented in figures 14-15. For our calcula- 
tions we assume that t in ~ 10 6 yr and note that results can be easily generalised for other 
values of t in . In figure 14 we show results of calculations of t ev for the two planetary masses 
(m p i = lMj, 5Mj), and two values of the angular velocity of the planet (Q r = 0, Vl crit ). The 
orbit of non-rotating lMj planet can be evolved by dynamical tides in a time < 10 10 yr for 
P bs < Pmax ~ 5days (see the lower solid curve in figure 14 ). The orbit of a lMj planet with 
Q r = Q crit has Pmax ~ 4:days. Orbits of non-rotating and rotating planets with m p \ = 5Mj 
have P m ax ~ 5 — Qdays. In that case the evolution is determined by dynamical tides in the 
star for P b s > 3days (both dashed curves almost coincide in figure 14 for P b s > 3days). 
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P obs (days) 



Figure 14. The time t ev (see equation (106) for definition) is shown as a function of P b s - The solid and dashed curves 
correspond to m p i = lMj, 5Mj, respectively. The lower and upper curves of the same type correspond to (Q r = 0, Q cr n). We 
choose cii n = lOAu for all curves. 




One can see from figure (14) that for a non-rotating planet and P b s ~ Iday, t ev is about 
a few years and is of order of the orbital period. Such a fast tidal evolution may possibly 
lead to disruption of the planet (see also Discussion). 
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t(10 6 yr) 

Figure 16. Same as figure 15, but for m p i = 5Mj. Note that all curves except the curve corresponding to U r = and 
Pobs = 3days are close to straight lines. Also, curves with the same period P j, s are close to each other. This is due to the fact 
that the orbital evolution is mainly determined by dynamic tides in the star for m p i = 5Mj . 

3.5 Onset of the stochastic instability 

Finally, let us discuss the condition leading to stochastic build-up of the mode energy. From 
the results of Section 2.9 it follows that the stochastic instability sets in when the orbital 
semi-major axis is sufficiently large: a > a st , where a st can be expressed in term of the 
parameter d crit (see equation (87)) using equation (84) 

a st = L - * C ; tEp \ F J 2/5 (— ) 3/6 fl* - 30.8(^ ooens )- 2 / 5 (^) 3 / 5 (^)^, 107 
QiTLU 00 (AE ns + AE ns *) m p i m pl Rj 

and we use the value a CT n = 0.67 on the right hand side. 

The result of calculation of a st for a IMj planet is shown in figure 17. One can see from 
this that the 'stochastic scale' is of the order of the scales of interest ~ 0.1 — WOAu, with 
larger values corresponding to the planet with Q r = Vt CTit . 

Accordingly, the stochastic build-up of mode energy may be difficult to achieve for planets 
rotating at the critical rate (see however Section 4). For non-rotating planets this problem 
is much less severe. For example, for the non- rotating planet with R p i ps 1.6Rj and P Q b s ~ 
Mays we have a st ~ 2Au (figure 17). Thus an orbit with 'initial' value of the semi-major axis 
~ lOAu can be evolved by stochastic tides in a time t ev ~ 10 5 yr (see figure 14) which is much 
smaller than the planet's cooling time ~ 10 6 yr. Therefore dynamical tides can significantly 
change the semi-major axes of non-rotating planets with parameters appropriate for observed 
systems. 

It is instructive to compare the analytic condition (107) with the results of numerical 
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P obs (days) 



Figure 17. The dependence of a s t (in Ait) on P b s (in days). The solid and dashed curves correspond to f! r = 0,Q cr n. Six 
curves of the same type correspond to different R p i/Rj = 1, 1.2, 1.4, 1.6, 1.8,2. The smaller values of a s t correspond to larger 
values of R pi . 

experiments (Mardling 1995 a, Mardling & Aarseth 2001). Madling (1995)a and Marldling & 
Aarseth (2001) define the boundary between stochastic and regular regimes of tidal evolution 
as a curve D^f n (e) on the plane (e,D min ), where e is the eccentricity of the orbit. The 
stochastic region is situated below this curve. For a n — 1.5 polytrope Mardling (1995)a 
obtains: D^Je « 0.9) « 3.75, D%g(e « 0.95) « 4.2, D^j*(e « 0.99) « 4.5, where the 
minimal separation or pericentre distance is expressed in units of stellar radii. Assuming 
that the mass ratio is equal to one and using the values lj o ~ 1-45 and Q m 0.5 (Lee & 
Ostriker 1986) we can obtain a similar curve D^ n (e) from our criterion (107). This curve 
is qualitatively similar to D£j*(e). We have D c ^ n (e w 0.9) w 3.1, £>^*(e w 0.95) w 3.4, 
Dmin( e ~ 0.99) ~ 4. Our curve lies below D^ n {e) on the plane (e, D min ) and for a specified 
value of e > 0.9 gives a value of -D m j n that is 10 — 20 percent smaller. A possible reason for 
this difference is the neglect of higher order terms in the tidal potential and the contribution 
of p-modes in our treatment of the problem which would enhance the tidal interaction and 
lead to larger values of D min . 

We do not discuss the condition for stochastic build-up of mode energy for planets with 
m p i = 5Mj. As was shown above, the semi-major axis of such a planet is mainly evolved as 
a result of dynamical tides in the star. These tides mainly excite / and low order g modes. 
These modes may interact effectively with high / g-modes due to resonant non-linear mode- 
mode interaction (Kumar & Goodman 1996). Kumar and Goodman estimated the energy 
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in tidally excited modes needed to be larger than ~ 10 36 ergs for the resonant parametric 
instability to occur. This is of the order of the energy obtained by the star after a fly-by 
of a 5Mj planet with minimal separation distances typical of our problem. The parametric 
instability may not only damp effectively the tidally exited modes but also lead to phase 
changes in these modes thus facilitating the stochastic build-up of mode energy in the star. 
However, this possibility needs further investigation. 



In this paper we have developed a theory of the tidal interaction between a central star and 
a massive convective planet moving on a highly eccentric orbit. We showed that it is possible 
to develop a self-consistent perturbative approach to the theory of disturbances induced in a 
slowly rotating planet by tidal interaction by regarding them as resulting from a sequence of 
impulsive interactions occurring at periastron (e.g.. Press & Teukolsky 1977). We obtained 
contributions to the energy and angular momentum exchanged during periastron passage 
from both quasi-static and dynamic tides that were derived from the same set of governing 
equations. Quasi-static tides require the existence of dissipative processes in order to produce 
net energy and angular momentum transfers from the orbit. However, dynamic tides may 
do this through the excitation of normal modes. 

We found that if the planet was initially non rotating, dynamic tides would transfer 
angular momentum to it until a 'critical' equilibrium prograde rotation rate, Vt crit , was 
attained. This is defined by the condition that dynamic tides acting on a planet rotating 
with the angular velocity fl cr it at periastron result in no net angular momentum transfer 
and hence do not change the rotation rate of the planet T . We found that the rotation rate 
of the planet which minimised the energy gained as a result of a periastron passage was also 
equal to Vt crit . We showed further that this rotation rate depends on the periastron distance 
between planet and star and that it can be larger than the similar equilibrium 'pseudo- 
synchronisation' rotation rate Q ps obtained through the action of quasi-static tides. In fact, 
for sufficiently small periastron distances, our analysis indicates that Q C rit can attain values 



of the order of the planet characteristic frequency = J „/' so that a planet rotating with 



Q r ~ £l C rit may be near to rotational break up. However, the analysis assumed slow rotation. 



4 DISCUSSION AND CONCLUSIONS 




^ Obviously, the value of Cl cr n would change 



in a more detailed treatment of the problem that takes account the contribution 



of other modes and the influence of quasi-static tides. 
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Thus this issue requires further investigation and consideration of the tidal interaction when 
the planet rotates rapidly. 

We examined the consequences of multiple periastron passages and energy exchanges, 
which, because they occur impulsively at periastron, could be modelled with a simple al- 
gebraic map. Multiple encounters can lead to a stochastic build-up of the energy contained 
in normal modes of oscillation (see Mardling 1995 a,b). We found a simple criterion for the 
instability of the dynamical system to such a stochastic build-up. 

Our results can be applied not only to recently discovered planetary systems containing 
a star and massive planet, but also in the much more general context of the tidal interaction 
of a solar type star with a fully convective object or the tidal interaction between two fully 
convective objects (such as e.g. low mass stars, brown dwarfs, etc.). 

We focused on the problem of the tidal evolution of a planet orbiting around a star on 
a highly eccentric orbit. Such a configuration might be produced as a result of gravitational 
interaction in a many planet system (e.g. Weidenschilling & Mazari 1996; Rasio & Ford 
1996; Papaloizou & Terquem 2001; Adams k Laughlin 2003). 

The minimum periastron distance of the orbit was assumed to be > 4 stellar radii 
corresponding to a final period P & s , after an (assumed) stage of tidal circularization at 
constant orbital angular momentum, exceeding about 3days. The initial semi-major axis a 
was assumed to be of the order of some 'typical' planet separation distance ranging from 
1 - lOOAu. 

We made some simple estimates of the possible tidal evolution of the orbit. We used a realistic 
model of the planet taking into account the fact that the radius and luminosity decrease 
with age as the planet cools. We first showed that quasi-static tides cannot account for a 
significant tidal evolution when the semi-major axes exceed ~ 0.1 Au. However, dynamical 
tides can in principle lead to a large decrease of the semi-major axis and circularize the orbit 
on a time-scale t ev ^ 10 9 — 10 10 yr. 

Note that for simplicity we did not take into account possible gravitational interaction of 
the planet with orbital parameters evolving due to tides with other planets in the planetary 
system. The possible characterization of this interaction could be decribed as follows. If the 
characteristic time of gravitational interaction t gr is smaller than t ev the orbital evolution 
would be mainly determined by the interaction with other planets. The interaction time t gr 
may increase significantly if the planet is scattered into a highly eccentric orbit by another 
planet that is itself ejected or scaterred to large distance. In addition t ev sharply decreases 
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with decrease of orbital semi-major axis. For a sufficiently small semi-major axis one would 
expect that t ev < t gr and the orbital evolution to be mainly determined by tides. Such a 
formulation of the problem is similar to the well-known loss cone problem in physics of stellar 
systems around super-massive black holes and is sufficiently complex to require analysis by 
numerical means. 

In a fully convective planet dynamic tides cause the excitation of mainly the / = 2 
fundamental mode of pulsation. This mode has a period that is typically much smaller than 
a characteristic time of periastron passage corresponding to a typical orbital period ~ P b s . 
The amount of energy and angular momentum gained at periastron due to dynamical tides 
decreases exponentially with P obs . Therefore, the evolution time t ev is very sensitive to the 
value of P obs . Orbits of planets with mass ~ lMj can be effectively circularized by the tides 
for final orbital periods P t, s ~ 3 — Adays corresponding to the shortest observed periods 
in extra-solar planetary systems. However, the contribution of dynamic tides to the orbital 
circularization rate does not seem to be effective for larger periods P b. s ^ 5days. These 
results appear to be consistent with the observational result that the orbits of extra-solar 
planets are nearly circular (eccentricity < 0.06) for periods P t, s ^ hdays. Dynamic tides 
raised in the star are not important for the orbital evolution of lMj planets. However, for 
planets with larger mass > 5Mj, the stellar dynamic tides are significant. 

There are at least three possible problems with the scenario of tidal circularization we 
have discussed here. Firstly, dynamical tides can only transfer energy from the orbit to 
the planet effectively only for sufficiently large orbital semi-major axes a > a st . For orbits 
with a < a st , stochastic instability leading to the energy transfer is ineffective so that the 
orbital evolution would instead be determined by the decay time of the fundamental mode 
of oscillation. This would lead to a very large circularization time-scale. This difficulty can 
possibly be overcome in the following way. Suppose that the planet starts its tidal evolution 
with a > a st . When a has decreased to a ~ a st , assuming little dissipation, a significant 
energy ~ GM j™vi would De contained in the oscillations of the planet. Eventually non linear 
effects should start to cause dissipation of this energy and a significant expansion of the 
radius which in turn produces a significant increase in the energy transfer from the orbit to 
the planet by tides. This could produce a form of tidal runaway (e.g.. Gu, Lin & Bodenheimer 
2003) until possibly a balance between dissipation of energy due to nonlinear effects and tidal 
injection is attained. In this situation, a condition of effective stochastic instability may be 
fulfilled all the way down to the final quasi-circular orbit. Self-consistent calculations taking 
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into account the non-linear effects of tidal heating on the structure of the planet are needed 
in order to evaluate this possibility II. 

On the other hand, if the dynamical tides operate effectively on scales <C a st there is 
enough energy exchanged to destroy the planet. Indeed if the planet semi-major axis is 
smaller than a* = ^-R v i ~ 7 • 10 12 (^)cm, the magnitude of the orbital energy becomes 
larger than the internal energy of the planet. Therefore, in order to settle onto a tight circular 
orbit around the star, short period planet must dissipate and radiate away an amount of 
energy which is much larger than its internal energy. This issue needs to be discussed in the 
context of a non-linear analysis of the orbital evolution. 

Finally, in our approximations, the influence of dynamical tides decreases exponentially 
with P obs and the tidal evolution as well as initial conditions of the problem are likely to be 
stochastic. Under such conditions there may exist other planets in the same system which 
have failed to circularize their orbits and have very high orbital eccentricities at the present 
time (the difference 1 — e ~ 10~ 3 — 1CT 4 ). This difficulty may appear less significant when one 
considers the possible contribution of other non-standard long period modes of oscillation. 
The energy transferred to these modes could lead to tidal evolution at larger separation 
distances (corresponding to P Q b s > 5). 

There are at least two possible candidates for such modes. Firstly, there is a low-frequency 
mode associated with the first order phase transition between molecular and metallic Hy- 
drogen (Vorontsov 1984, Vorontsov, Gudkova and Zharkov 1989). This mode has a period 
approximately twice as long as the period of the fundamental mode. However this mode is 
probably not important for our problem. As we have mentioned before the phase transition is 
absent in hot young planets. Also, this mode disappears when the time scale of the transition 
of one phase to another is smaller than the period of the mode (Vorontsov 1984). Another 
very interesting possibility is to consider the inertial modes associated with rotation of the 
planet (e.g. Papaloizou & Pringle 1978). The modes have very low frequency and therefore 
their excitation does not decrease exponentially with distance. The density perturbation 
associated with these modes is small being of order fi 2 .. However, the contribution of these 
modes to the energy gained from the orbit may be significant for large values of P b s where 

II Tidal heating in a similar context has been discussed by Podsiadlowski (1996) for stellar systems and more recently by Gu, 
Lin & Bodcnhcimcr (2003) for short period planets heated by the action of quasi-static tides. 



© 1994 RAS, MNRAS 000, 1-42 



40 P. B. Ivanov, J. C. B. Papaloizou 

the contribution of the fundamental mode is exponentially small. We will to consider the 
excitation of these modes by dynamical tides in our future work. 
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APPENDIX A: DERIVATION OF THE EXPLICIT EXPRESSION FOR 
THE DAMPING RATE 

In order to obtain equation (20) from equation (15), we need to evaluate the following 
expression 

a = \ J &lo*o$, Al 

where o*® is related to the displacement through equation (16) and integration is performed 
over the solid angle Q. We assume an isentropic condition leading to a circulation free 
displacement such that V x £ = is fulfilled. We may then introduce a displacement potential 
$ such that £ = V$, and rewrite (Al) in terms of it so obtaining 

a = 2j dn($* a ^, af3 - W A$). A2 

We use the identity 

®: a ^,«p = - ($; a A$ Q + $ q a$;j), 

to bring (A2) in the form 

a = J ^(A(<T Q $ Q ) - (<T a A$ Q + $ q a$;j _ jWA$) 

= J dQ(A(£ ■ O - C ■ V(V • - £ • V(V • D - |(V • 0(V • £*))■ ^3 

Now we substitute equation (12) into (A3), use the relation (13) and the known properties 
of the vector functions to obtain 

a = 2{ l^ R ~V + Y7 )2 + ~ T^ 2 + ^ - 2 ^ + L ^ R - ^ M 

where L 2 = 1(1 + 1). Note that (A4) cannot be used for I = 0. Multiplying equation (A4) on 
(r 2 pz/), setting / = 2 (L 2 = 6) and integrating the result over r we get equation (20). 
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